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INTEGRATED ORBIT AND ATTITUDE CONTROL FOR A 
NANOSATELLITE WITH POWER CONSTRAINTS 

Bo J. Naasz,* Matthew M. Berry,* Hye- Young Kim,* 
and Christopher D. Hall§ 


Small satellites tend to be power-limited, so that actuators used to 
control the orbit and attitude must compete with each other as well 
as with other subsystems for limited electrical power. The Virginia 
Tech nanosatellite project, HokieSat, must use its limited power re- 
sources to operate pulsed-plasma thrusters for orbit control and mag- 
netic torque coils for attitude control, while also providing power to 
a GPS receiver, a crosslink transceiver, and other subsystems. The 
orbit and attitude control strategies were developed independently. 
The attitude control system is based on an application of LQR to an 
averaged system of equations, whereas the orbit control is based on 
orbit element feedback. In this paper we describe the strategy for 
integrating these two control systems and present simulation results 
to verify the strategy. 


INTRODUCTION 

The concept of orbit and attitude coupling has received little attention in the literature. 
The time scales of these two systems usually differ enough that the orbit can be considered 
as given when considering attitude dynamics. In some cases, however, this assumption is 
unwise, and consideration must be given to the coupling of the translational and rotational 
systems. Integration of spacecraft orbit and attitude subsystems is motivated by a variety of 
concepts including natural system dynamics, and guidance, navigation and control (GNC) 
system issues such as shared resources, subsystem inter-dependencies, and actuator-induced 
disturbances. 

Natural coupling of rotation and translation is apparent in such concepts as variation 
of atmospheric drag and solar radiation pressure as a function of attitude-dependent cross- 
sectional areas, variation of magnetic and gravity gradient torques as a function of both 
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attitude and orbit position, and dynamical coupling of orbit and attitude dynamics (a weak 
effect which is most often negligible). 1,2 

There are numerous examples of shared resources, including shared actuators (thrusters 
for orbit and attitude control, as well as for momentum dumping maneuvers), shared sen- 
sors (optical systems for attitude determination and celestial navigation), and other shared 
resources such as power and spacecraft consumables. Subsystem inter-dependencies include, 
for example, the dependence of the orbit control system (OCS) on the attitude control sys- 
tem (ACS) for actuator pointing. This concept is especially important for spacecraft em- 
ploying atmospheric drag or solar radiation pressure control techniques. Actuator-induced 
disturbances include effects such as torques and forces induced by the use of non-coupled 
or misaligned thrusters to perform rotational or translational maneuvers. This concept also 
introduces the need for cooperation between attitude and orbit control and determination 
subsystems to support accurate state estimation. 

In this paper, we explore some of these orbit-attitude integration issues as applicable to 
the HokieSat spacecraft. 3 The spacecraft is one of three small spacecraft in the Ionospheric 
Observation Nanosatellite Formation (ION-F). 4 The overall configuration and structural 
design of the satellite is described in Refs. 5 and 6. The attitude control system, including 
the hardware configuration and the algorithm development, is described in Refs. 7, 8, and 9. 
The orbit control system is described in Refs. 10, 11, and 12, and related optimal control 
concepts are described in Refs. 13 and 14. 

MISSION DESCRIPTION AND HARDWARE CONSTRAINTS 

HokieSat is a small spacecraft being built by students at Virginia Tech. 3 The spacecraft, 
formally known as the Virginia Tech Ionospheric Scintillation Measurement Mission, is part 
of the Ionospheric Observation Nanosatellite Formation (ION-F) project. 4 The project 
involves spacecraft built by three schools: Virginia Tech, Utah State University, and Cornell 
University. The three spacecraft are similar in design and will perform formation flying 
demonstrations and make ionospheric measurements while in an orbit approximately the 
same as that of the International Space Station ( i « 52°, a « 6770 km). 

HokieSat is hexagon-shaped, about 0.5 m in major diameter and about 0.3 m tall, with 
a mass of about 20 kg. The spacecraft uses Pulsed Plasma Thrusters (PPTs) to maintain 
its position in the formation. 15,16 There are two pairs of PPTs on HokieSat; their position 
on the hexagonal cross-section is shown in Fig. 1. Thrusters T2 and T3 provide translation 
control, and Thrusters T\ and T4 can provide yaw steering, or can provide translation 
control. Each thruster can be fired individually However because they share a capacitor, 
thrusters T\ and T2 or thrusters T3 and T4 cannot be fired simultaneously. Thrusters T2 
and T3 can be fired simultaneously, as can thrusters T\ and T4. Each thruster provides an 
impulse-bit of 56 /xN-s and fires at a rate of 1 Hz. For translation control thrusters T 2 and 
T3 are fired together providing an impulse-bit of 112 ^xN-s. All four thrusters are positioned 
slightly above the center of mass, and therefore exert a small pitch or roll torque on the 
spacecraft. There are no thrusters in the zenith-nadir directions, and the communication 
system requires that the spacecraft remain nadir-pointing; thus there is no way to thrust in 
the radial direction. 

Attitude is controlled by three orthogonal magnetic torque coils. 7,8 Attitude control is 
achieved by forcing a current through the torque coils, thereby generating a net magnetic 
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Figure 1: HokieSat pulsed plasma thruster (PPT) layout viewed from zenith direction 

moment which interacts with the Earth’s magnetic field and creates a torque. Due to 
magnetic field interactions between the coils and the PPTs, the two actuator systems cannot 
be used simultaneously, and any attitude or orbit control must be performed in a piecewise 
fashion. Power limitations place an additional constraint on the control system. When the 
spacecraft is in eclipse, the power subsystem can provide only enough power to operate 
vital spacecraft functions; therefore, attitude control can be accomplished during eclipse, 
whereas orbit control can only be accomplished when the spacecraft is in sunlight. 

HokieSat ’s nominal flight orientation is nadir pointing, with thrusters T2 and T3 pointing 
opposite the velocity vector (see Figure 1). The most obvious control mode employs torque 
coils to slew the spacecraft as required to point thrusters T2 and T3 in the direction dictated 
by the orbit control law. This mode is termed the “normal” control mode. A less obvious 
control strategy is to turn the spacecraft sideways, so that thrusters T2 and T3 point in 
the negative orbit normal direction, with thruster T\ and T4 aligned with the velocity and 
negative velocity directions. This mode is termed the “sideways” control mode. 

The normal control mode benefits from two parallel thrusters (and thus doubled AU), 
and no thrust torque disturbance about the yaw axis. The primary disadvantage is the 
inherent requirement that the attitude control system perform 180° slew maneuvers when- 
ever the requested thrust vector changes directions. We can modify the orbit control law 
to discourage this thrust vector flipping by recognizing that all elements except the semi- 
major axis (and thus the true anomaly) can be controlled without thrust vector flipping. 
For example, positive changes in the inclination are created by thrusting in the orbit normal 
direction when the spacecraft is above the equator. If negative changes are required, the 
control law can simply wait for the spacecraft to move below the equator, where thrust in 
the orbit normal direction results in negative inclination change. This observation is use- 
ful. Unfortunately, there is no way to force positive and negative semi-major axis changes 
without applying thrust in opposing directions (the positive and negative velocity vector 
directions), and large slew maneuvers are unavoidable for the normal control mode. 

The sideways control mode benefits from the fact that thrust is readily available in the 
positive and negative velocity directions, as well as one out-of-plane direction. The major 
disadvantage of the sideways control mode is the torque induced by orbit control thrusting. 
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Thrusters T\ and T4 are not aligned with the center of mass, and thus create yaw torque 
about the zenith-nadir directions. This disturbance torque has the same order of magnitude 
as the torque available from the torque coils. In this control mode, the out-of-plane elements 
are controlled as described above, with thrust only in the positive orbit normal direction. 
The sideways control mode requires no large slew maneuvers, but does require attitude 
control to compensate for the thruster disturbance torque. 

Reference 10 described various linear and nonlinear control laws for HokieSat to main- 
tain the formation. These orbit control laws include the non-radial thrust constraint, and 
the constraint on the thrust magnitude. The control laws assume that the necessary atti- 
tude can be achieved instantaneously, and that there is no power constraint. Reference 7 
described an LQR-based attitude control for the HokieSat mission. The control law was de- 
veloped to maintain nadir-pointing, but not to perform yaw-steering for thruster pointing. 
The control law did not enforce any power constraints on torque coil operations. 

The coupling of these power-constrained orbit and attitude control laws is an interesting 
nonlinear control problem. Interruption of the orbital control may drive the system unstable. 
In this paper, we implement integrated nonlinear orbit and attitude control laws and test 
their performance in a realistic formation simulation. 

EQUATIONS OF MOTION 

In this section we describe the system dynamics models used in verification of the control 
strategies. 

Keplerian Orbital Dynamics 


In cartesian form, the equations of motion for a point-mass satellite take the form 


r = — - 


rr + a p 


(i) 


where r is the position vector of the satellite measured from the center of mass of the 
primary body, ji is the gravitational parameter of that body, and a p includes perturbation 
accelerations caused by oblateness effects of the body, atmospheric drag, third body effects 
and so on. If a p = 0 then the equations of motion reduce to the ideal Keplerian situation 
where the body is a perfect sphere and all other disturbances are zero. 

We use the Earth-centered inertial frame, with the i3-axis aligned with the north 
pole, the ii-axis aligned with the equinox direction through the equatorial plane, and the 12- 
axis completing the orthonormal set. We also refer to the spacecraft orbital frame, T 0 , which 
is centered on the spacecraft, with the 03-axis aligned with the vector from the spacecraft to 
the center of the Earth, the 02-axis aligned with the negative orbit normal direction, and the 
61-axis aligned with the velocity direction for circular orbits. The spacecraft body frame, 
Tb, is centered on and fixed to the spacecraft. The bi-axis is in the normal mode thrust 
direction (the direction opposite thrusters T2 and T3), the 63 direction points through the 
nadir panel, and the b2 direction completes the right-hand set. 

We define the classical orbital elements as semi-major axis, a, eccentricity, e, inclination, 
i, right ascension of the ascending node, fi, argument of periapse, cj, and mean anomaly, 
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M. The first five orbital elements describe the size, shape, and orientation of the orbit 
in inertial space. The final element describes the spacecraft’s angular position within the 
orbit, referenced from periapse, and has constant angular rate given by the mean motion, n. 
The equations of motion of a controlled spacecraft in terms of the classical orbital element 
set are given by Gauss’s form of Lagrange’s planetary equations: 17 
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where u r , and u & are the radial, transverse, and orbit normal control acceleration 
components, v is the true anomaly, # = cj 4* v is the argument of latitude, and 6, h, and p 
are the semi-minor axis, the angular momentum, and the semi-latus rectum, respectively. 


Attitude Dynamics 


We use Euler parameters, or quaternions, to describe the spacecraft attitude. The quater- 
nion, q is defined in terms of the Euler axis, a, and the Euler angle, 


. $ 
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( 9 ) 

(10) 


The time rate of change of the quaternion is written in terms of the angular velocity, cj, as 


q =2 


q x + 94 1 
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( 11 ) 


where q x is the skew symmetric matrix: 
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( 12 ) 


The rotational equations of motion for a rigid body, Euler’s equations, are 

I Cj =■ -u x Iuj + g e 


( 13 ) 
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where I is the moment of inertia matrix, and g e is the external torque on the body. 

For an orbiting spacecraft, external torques include the control torque, g c , the gravity 
gradient torque, g gg , and other lesser torques from atmospheric drag, solar radiation pres- 
sure, and other environmental effects. For the simulations presented here, we consider only 
control torques (both intended and unintended) and gravity gradient torques. The gravity 
gradient torque is given by 

g 99 = 3Wc°3 X I°3 (14) 

where uj c is the mean motion (also denoted n). 

ATTITUDE CONTROL 


The attitude control torque is effected by controlling the current in three mutually orthogo- 
nal magnetic torque coils, thereby generating a magnetic moment, M. The resulting torque 
is g — M x B, where B is the magnetic field vector. The B field varies with both orbital 
and diurnal periods, and it is clearly not possible to cause a torque about the B direction. 
Thus, magnetic attitude control is inherently underactuated and involves time-varying pa- 
rameters. The approach taken in Refs. 7 and 8 was based on linearizing the attitude motion 
about the nominal nadir attitude, averaging the B field, and applying LQR to the resulting 
linearized equations of motion. Subsequent Floquet analysis and nonlinear simulations were 
used to verify performance of the resulting controller. For further details, see Refs. 7 and 8. 

ORBIT CONTROL 


The orbital equations of motion (2 - 7) can be written as: 

de = f(oe) + G(oe)u (15) 

where G(oe) is the input matrix, whose elements are evident from Eqs. (2-7). We separate 
these equations of motion into two systems: one system consisting of the first five elements, 
which define the shape, size, and orientation of the orbit, and the second system consisting 
of the in-plane angular position of the spacecraft. 

For the first system, we define the elemental error 77 as 
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where * designates target elements. The equations of motion for this system are: 


(16) 


r?= Gu 


(17) 


where u is the vector control acceleration in radial, transverse, and orbit normal directions, 
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and the input matrix, G comprises the first five rows of G(oe): 
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where the sine and cosine functions are abbreviated as S(-) and C(-). A candidate Lyapunov 
function for this system, weighted by positive gains, Kj. is given by 


V{rj) = i [K a 8a 2 + K e Se 2 + Krfi 2 + K^SQ 2 + KJlo 2 ] 


The time derivative of V ( 77 ) is given by 

V(ri) = V v ri = V f ,G(ri)u 

where 

= [ K a Sa K e 8e KiSi K Q 5tl K^Slj 
Choosing the control, u, as 


u 


= -g^vi 


results in a negative semi-definite time derivative of V: 

V(fl) — - Vt?GG t V^ 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


To prove asymptotic stability, we apply LaSalle’s invariance theorem. The time deriva- 
tive of the Lyapunov function is always zero when 77 = 0 , and could be zero when the 
trigonometric functions of 6 are zero, which occurs when 9 = kit/ 2, where k is an integer. 
The set where 9 = kn/2 is not an invariant set, because 9 is time- varying and therefore 
trajectories that start in the set do not remain in the set. Therefore, the largest invariant 
set where V = 0 is the set 77 = 0 . By Lyapunov’s direct method, and LaSalle’s invariance 
principle, the system is asymptotically stable under the choice of control given in Eq. (24). 
Furthermore, the system is globally asymptotically stable because V — > 00 as ( 77 ) —> 00 . 
Thus, the desired feedback control is given by 
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The angle errors SQ, Su), and SM are measured from parameters defined between 0 and 
27 t. To ensure proper error feedback, these angular errors should be “short-way” angle 
measurements, defined between -n and n. 
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Mean Motion Control 


To control the in- plane angular motion of the spacecraft, we take advantage of a useful 
natural component of the nonlinear dynamics: the differential mean motion, Sn. Observe 
that an effective method for correcting in-track errors is to force a semi-major axis error, 
changing the mean motion such that the in-plane angle is corrected. For an uncontrolled 
spacecraft in two-body motion, the orbital element dynamics simplify to 


dM 

dt 


= n 


(25) 


and thus the relative dynamics can be written as 


d f ..^ dM dM* , . 

M {SM) = ^-nr = n ^ n = Sn 


where 



To drive the mean anomaly error, SM, to zero, we use 


SM = ~K n SM 


(26) 

(27) 

(28) 


where K n is a positive gain. Combining Eqs. (27) and (28), and using canonical units 
((i = 1), we obtain 



Solving for a, we define a new target semi-major axis, a**, which forces the mean anomaly 
error to zero: 

/ 1 \-2/3 

a” = (- K n SM + (30) 

Notice that as SM goes to zero, the mean motion control target semi-major axis, a**, 
approaches the original target value, a*. In application, we replace the mean anomaly 
error, SM, in Eq. (30) with argument of latitude error, SO , so that the mean motion control 
properly positions the spacecraft within the orbital plane, even in the presence of argument 
of periapse error. 

Using the orbit control law defined in Eq. (24), with the target semi-major axis a ** de- 
fined in Eq. (30), we can control the full, nonlinear motion of an orbiting spacecraft. Since 
most spacecraft have fixed magnitude thrusters, we must develop some way of discretizing 
the control acceleration requested by the control laws. For example, if the orbital element 
feedback control law requests acceleration of magnitude A in the a direction, we must deter- 
mine whether or not to thrust, and if so, in what direction. This decision is straightforward 
for this elemental control law, as we have full scaling freedom in the choice of gains. We 
simply calculate a desired thrust direction and magnitude from the control law, and thrust 
if and only if the desired thrust magnitude is greater than the available thrust magnitude. 
Mathematically, this thruster on/off logic requests thrust in the a direction if and only if 
A > T/m where T is the available thrust magnitude, and m is the spacecraft mass. 

In the case of the HokieSat mission, we not only scale the thrust magnitude, but 
also constrain the thrust direction. HokieSat ’s orbit control consists of four pulsed plasma 
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thrusters (PPTs) aligned in the local horizontal frame, with no thrust available in the radial 
direction. The classical orbital element Lyapunov control law can be modified to exclude 
radial thrusting by setting the terms in the first column of the G matrix in the control 
definition, Eq. (24), to zero, thus requesting only transverse and orbit normal thrust. By 
LaSalle’s invariance principle, this modification does not affect the stability of the system, 
as the maximum invarient set for which V = 0 remains 77 — 0. 


Orbit Control Gain Selection 

Nonlinear feedback control gain selection techniques are described in Ref. 11. Following 
these guidelines, we choose the following positive feedback control gains to 1) guarantee 
asymptotic convergence to an error envelope, 2) appropriately weight the elemental errors, 
and 3) eliminate the problem of chattering by properly defining the size of the error envelope. 
These gains are calculated as functions of the target orbit elements in canonical units (/z = 1, 
DU — R e ) as follows: 
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We modify these feedback control gains by trigonometric functions of the true anomaly, v , 


as suggested by Schaub and Alfriend: 18 


K a = K a sin N v 

(36) 

K e = K e cos N v 

(37) 

ki = Ki cos N 6 

(38) 

Kq = K Q sm N 0 

(39) 

Ku = K u sin N v 

(40) 


where N is some positive integer. Schaub et ai 18 suggest a slight variation of this form 
(with K a multiplied by cos N z^/2), to encourage control of the elements when they are the 
most controllable. For example, semi-major axis error would be emphasized at periapse, 
and mostly ignored elsewhere. The form given by Eq. (36) is more appropriate for HokieSat 
because it encourages control of Sa not when a is most controllable, but when semi-major 
axis control will have the least impact on the eccentricity error. By not creating error in 
one element in the process of correcting error in another, 19 we can significantly reduce the 
fuel required to perform a given maneuver. 
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INTEGRATED ATTITUDE AND ORBIT CONTROL SIMULATION 


To test the control laws, we apply them in a simulation of HokieSat’s flight software. The 
test case is a rendezvous between a controlled spacecraft and an uncontrolled target orbit. 
The simulation is written in C++ and performed on a Linux machine. The simulation 
emulates the software environment that will run on HokieSat’s onboard computer. 

HokieSat will be controlled by an onboard computer running the VxWorks real-time 
operating system. The software architecture uses several processes running simultaneously. 
These processes can set global variables that other processes can read, and the processes can 
also send messages to one another. Each process controls a specific spacecraft subsystem; 
e.g., there is a process that controls downlink communications to the ground station. 

In these simulations, we model the software architecture that will be used onboard 
HokieSat. The simulations are performed in the Linux operating system, using a C++ li- 
brary that simulates the VxWorks operating system. The simulation involves four processes, 
which interact in the same way that the flight processes interact. The four processes are 
Orbit Determination, Orbit Control, Attitude Determination, and Attitude Control. The 
intention is for each of the controls processes to resemble the flight code as much as possible. 
The attitude and orbit determination processes are simulated by numerical integrators, in 
place of the sensors that will be on HokieSat. 

Orbit Determination Process 

At each time step, the orbit determination process propagates the orbital position and ve- 
locity from the previous time to the current time. The propagation is performed using a 
4th-order Runge-Kutta integrator. The only forces included are the two-body force and the 
applied thrust. Both orbits are propagated at each time step; the target orbit is integrated 
without any thrust, and the controlled orbit is integrated using the thrust direction calcu- 
lated by the orbit control process in the previous time step. The magnitude of the thrust 
is 56 per thruster used, or zero if no thrust is applied. We assume a 20-kg spacecraft. 

Attitude Determination Process 

The attitude determination process propagates the attitude of the controlled spacecraft at 
each time step. The equations of motion for the attitude are integrated with a 4th-order 
Runge-Kutta integrator to find the current quaternion relating body and orbital frames, 
and the angular velocity, at the current time. Thrust and gravity gradient torques are 
considered, as well as the magnetic control torque calculated by the attitude control process 
at the previous time step. 

Orbit Control Process 

The orbital control process calculates the control thrust at each time step. It reads the 
current and desired (target) orbital positions and velocities, and determines the desired 
thrust direction based on the control law. Because HokieSat must remain nadir pointing, 
any radial component of the thrust direction determined by the control law is ignored and 
set to zero. If the control law determines a thrust magnitude that is less than the amount of 
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thrust that HokieSat can provide, then no thrust is applied. If the control law determines a 
thrust that is greater than the amount of thrust that HokieSat can provide, then the thrust 
magnitude is set equal to HokieSat’s thrust. Based on the desired thrust direction, a desired 
attitude is determined in a quaternion form. This quaternion is compared to the current 
attitude quaternion calculated by the attitude determination process. If the attitude error 
is within 5°, and if the angular velocities are below 0.001 rad/sec, thrust is applied. When 
the thrust is applied, the thrust direction in the inertial frame is calculated by converting 
the thrust direction from the body frame to the inertial frame. This inertial thrust direction 
is then sent to the orbit determination task to use when propagating to the next time step. 
If the angle between the desired attitude and the current attitude is above the tolerance, or 
the angular velocities are above the tolerance, no thrust is performed, and the orbit control 
process sets the desired attitude for the attitude control process to read. 

Attitude Control Process 


At each time step, the attitude control process reads the current and desired attitudes and 
angular velocities, and calculates a torque based on the control law. If the calculated torque 
magnitude is above the maximum torque that HokieSat can provide, the torque is set equal 
to the maximum torque. If the orbital control process has not commanded a thrust, then 
the calculated torque is sent to the attitude determination process to use when propagating 
to the next time step. If the orbit control process has commanded a thrust, then the control 
torque is set to zero. 

SIMULATION RESULTS 


The initial conditions of the target orbit are defined by a 6770 km, nearly circular orbit with 
inclination of 52°. The initial conditions of the maneuvering spacecraft are dependent on 
the test case, and are either a 700 m position separation from the target, a 1 m/s velocity 
error from the target, or a combination of the two. In the simulations the initial attitude 
is chosen so that there is no initial attitude error. The initial angular velocity between 
the orbital and body frames is set to zero. The moment of inertia matrix, I, used in the 
simulation, is based on a mass properties study performed at Wallops Flight Facility, 6 and 
is: 


1 = 


0.518 

0.0046 

0.0 


0.0046 

0.4898 
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0.0 

0.0 

0.4530 


kg m 2 


(41) 


We investigate several possibilities using these simulations. There are five types of cases 
to consider for simulations, each with varying degrees of fidelity as regards the actual system. 
First, the attitude control torque can either be the ideal torque or it can be constrained to 
the actual magnetic torque coil system that will be used in the satellite. The latter is of 
course the most important; however, simulations using the ideal torque aid in understanding 
of the overall system performance. Second, the orbit control strategy can use the normal 
mode or sideways mode. Third, the spacecraft can execute required controls throughout 
the orbit or the effects of eclipse can be included. For the eclipse mode, we assume as a 
worst case that eclipse lasts 45 minutes. Fourth, the initial position error (taken as 700 m) 
can be either in-track, cross-track, or radial, or a combination of the three. Fifth, the initial 
velocity error (taken as 1 m/s) can be in any direction. 
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Figure 2 shows the results of a leader-follower maneuver performed in the normal flight 
mode, with magnetic torque attitude control, and no eclipse eclipse constraint. The initial 
state of the maneuvering spacecraft in this simulation is defined by a 1 x 10“ 4 radian initial 
mean anomaly offset from the target orbit. Position error is reduced to < 10 m in less 
than one day. The convergence rate of position error is slow because the orbit control 
law requires frequent 180° attitude maneuvers but the available magnetic attitude control 
torque is small. 


Position Error - S2001 0 x 1 0" 4 A PP lied Thrust - S2001 0 




Angle Error - S2001 0 x 1 0‘ 6 A PP |ied Torque - S2001 0 




Figure 2: Simulation results for a 700-m leader-follower maneuver in the normal flight mode. 

Figure 3 shows the results of a same-ground-track maneuver performed in the normal 
flight mode with magnetic attitude control, zero velocity error, and no eclipse. The initial 
state of the maneuvering spacecraft in this simulation is defined by a 1 x 10“ 4 radian initial 
mean anomaly offset, and a — 1 x 10“ 5 radian initial right ascension offset from the target 
orbit. In this simulation, the control fails to complete the maneuver. This is due to the 
rapidly changing thrust direction, and the magnetic attitude control subsystem’s inability 
to promptly point the thrusters in the required direction. 
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Figure 3: Simulation results for a 700-m same ground track maneuver in the normal flight 
mode. 

This problem due to the attitude control by the magnetic torque coil can be overcome by 
the sideways flight mode. Figure 4 shows a 700-m same ground track maneuver, respectively, 
in the sideways flight mode with the magnetic torque attitude control. The convergence 
rate of position error is much faster than the previous case because this mode does not 
require a large angle attitude maneuver. 

Considering the effects of eclipse has further impact on performance. In Figure 5, the 
eclipse lasts for a half orbit period per each orbit. This trajectory uses sideways mode, 
in-track position error, magnetic attitude control, and no velocity error. The convergence 
rate is slightly slower than the no eclipse case but not by much. Attitude control during 
the eclipse allows the orbit control to perform effectively immediately after eclipse. 

Figure 6 shows simulations results with initial velocity errors. When the initial velocity 
error is large enough, it is difficult to overcome these errors using a limited power thruster, 
especially if the velocity error is in the orbital velocity direction (Figure 6(a)). It takes 
hundreds of orbits to be able to observe convergence. When considering eclipse, it takes 
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Figure 4: Simulation results for a 700-m same ground track maneuver in the sideways flight 
mode. 

twice as long to achieve convergence. This result implies a significant emphasis should be 
placed on insuring that initial velocity errors due to spacecraft separation are minimized. 

CONCLUSIONS 

Attitude and orbit control systems are frequently considered to be independent. However, in 
small spacecraft with limited power, the coupling can be significantly greater if control actu- 
ators cannot be used simultaneously. Further performance degradation results if actuators 
cannot be used during eclipse. Independently developed controllers can be successfully inte- 
grated to achieve reasonable performance, as illustrated by the simulation results presented 
in this paper. The HokieSat mission has limited power, limited thrust magnitude, fixed 
thrust direction, no thrust during eclipse, and limited attitude control torque. A successful 
control strategy uses orbit element feedback for orbit control, and a modified LQR-based 
control for attitude control. To avoid the problems of frequent direction changes in the re- 
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Figure 5: Simulation results for a 700-m leader-follower maneuver in the sideways flight 
mode with eclipse. 

quired thrust, the spacecraft flies “sideways” so that thrust can be applied in the ± velocity 
direction. Effects of initial velocity and position errors can be substantial. 
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Figure 6: Simulation results for initial velocity errors (lm/sec) in the sideways flight mode 
(a) Oi-axis direction error - S21002, (b) random direction velocity error - S21003, (c) position 
and random direction velocity error - S21013, (d) oi-axis direction error with eclipse - 
S21102. 
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